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Abstract.  Large  eddy  simulation  (LES)  with  a flarnelet-based  chemistry 
model  has  been  evaluated  through  a priori  and  a posteriori  tests  in  both  de- 
caying homogeneous  turbulence  and  spatially-developing  mixing  layers.  The 
present  flamelet-LES  approach  involves  a series  of  models,  among  which 
are  the  models  for  the  PDF  of  mixture- fraction  (/  ( z ))  and  the  conditional 
filtered  dissipation  rate  (Xst)-  We  test  three  models  for  f (z):  a)  6— PDF 
(/  (z)  = 6 (z)),  b)  /3-PDF,  c)  Gaussian  PDF;  and  two  models  for  Xst • The 
Gaussian  PDF  model  consistently  performs  almost  as  well  as  the  /3— PDF, 
and  may  provide  an  alternative  approach  to  the  (3— PDF  method.  This 
point  is  important  to  the  LES  simulations  since  the  calculation  of  Gaussian 
function  is  much  cheaper  than  that  of  /3-function.  Furthermore,  it  may  be 
said  from  the  error  analyses  that  the  counterflow  model  for  Xst  has  average 
performance  in  its  present  form  and  that  a phenomenological  model  in  the 
form  of  Eq.  (6)  exists  with  an  optimum  value  of  a parameter  C3  to  produce 
a better  result.  The  a posteriori  analysis  shows  a satisfactory  performance 
of  the  flamelet  model  within  the  context  presented  in  this  paper. 


1.  Introduction 

Large  eddy  simulation  (LES)  constitutes  an  attractive  approach  for  nu- 
merical simulation  of  turbulent  reacting  flows.  The  basic  idea  of  LES  is  to 
calculate  the  large-scale  energy-containing  eddies  and  use  a subgrid  model 
for  the  small  scales.  The  large-scale  structures  resolved  by  LES  are  effec- 
tive entrainers  and  play  a role  in  bringing  various  reactive  gas  pockets  into 
contact  before  the  reactants  are  mixed  by  molecular  diffusion,  prior  to  reac- 
tion. Hence,  the  entrainment  rates  induced  by  the  large-scale  structures  will 
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determine  the  overall  reaction  rates  in  a turbulent  reacting  flow,  and  are 
sometimes  crucial  to  the  understanding  of  flame  behavior,  especially  when 
combustion  instabilities  occur1.  LES  is  thus  regarded  as  a favorable  tool 
in  combustion  applications,  better  than  the  traditional  Reynolds- Averaged 
Navier-Stokes  (RANS)  turbulent  models  which  are  limited  to  the  descrip- 
tion of  the  mean  flow  field. 

However,  relatively  few  studies  have  addressed  LES  of  reacting  flows. 
Modeling  subgrid  combustion  activities  encounters  great  challenges.  Un- 
like the  aerodynamic  problem,  the  use  of  the  similarity  assumption  for 
the  small-scale  mixing  and  dissipation  processes  in  reactions  leads  to  un- 
resolved terms  which  are  related  to  the  heat  release.  In  addition,  a more 
serious  problem  for  LES  combustion  models  is  that  the  chemical  reactions 
almost  always  take  place  within  the  unresolved  subgrid  scales.  For  example, 
an  approximate  model  neglecting  subgrid  scale  contributions,  i.e.,  writing 
the  reaction  rates  as  an  Arrhenius  law  in  terms  of  the  filtered  quantities, 
significantly  misrepresents  the  combustion  dynamics2.  The  modeling  tasks 
are  then  as  challenging  as  in  RANS  applications. 

The  current  practice  in  modeling  the  subgrid  combustion  activities  gen- 
erally follows  the  concepts  and  the  techniques  once  developed  for  RANS 
applications.  They  include  the  direct  methods,  such  as  the  extended  version 
of  the  Eddy-Break-Up  model2,  the  Linear- Eddy  model3,  the  Transported 
Probability  Density  Function  (TPDF)  method4  and  the  Conserved  Scalar 
method,  such  as  the  flamelet  approach5.  Some  new  models  based  on  the 
similarity  concept  have  also  been  suggested6.  Until  recently,  the  focus  was 
on  a priori  testing  of  the  applicability  of  combustion  models  in  LES7;  the 
evaluation  of  the  models  by  a posteriori  testing  is  not  as  common,  except 
in  the  studies  of  Menon  and  his  co-workers8  and  Pitsch  and  Steiner9. 

The  basic  procedure  for  the  flamelet-LES  model  used  in  the  present 
paper  is  contained  in  Cook  and  Riley7.  However,  the  present  paper  uses 
this  technique  within  the  framework  of  a generalized  curvilinear  coordinate 
system  to  permit  the  calculation  of  turbulent  combustion  in  realistic  sys- 
tems which  usually  have  complex  geometries.  The  main  contribution  of  the 
present  paper  is  in  the  investigation  of  the  various  model  assumptions  used 
for  the  calculation  of  the  mixture-fraction  dissipation  rate  and  its  PDF. 
Both  the  a priori  and  a posteriori  testing  through  direct  numerical  simu- 
lation (DNS)  and  LES  of  turbulent  non-premixed  flames  are  reported. 

2.  Subgrid  Flamelet  Models 

The  generalized  curvilinear  coordinate  formulation  for  the  large-eddy  sim- 

ulation follows  Jordan11,  where  the  full  Navier-Stokes  equations  are  trans- 

formed prior  to  filtering.  The  filtered  reaction  rate  is  approximated  by  the 
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laminar  flamelet  model.  That  is,  in  the  unresolved  subgrid  scales,  combus- 
tion takes  place  in  a thin  layer  in  the  vicinity  of  the  surface  of  stoichiomet- 
ric mixture  fraction  where  its  local  gradient  is  sufficiently  high.  Therefore, 
the  combustion  in  a turbulent  flow  is  represented  by  a statistical  ensemble 
of  such  laminar  flamelets.  At  each  grid  point,  the  filtered  reaction  rate  is 
therefore  modeled  as 

5i(x)  = J f u>i{Da,Ze,Ce,Po,YiuYi2,TuT2]Z,Xs)  ■ f {Z,Xs;*)dZdxs- 

1(a) 

The  mass  fraction  is  also  modeled  in  a similar  fashion: 

Yi  (x)  = J I Yi(Da,Ze,Ce,P0,Yii,Yi2,Tx,T2-,Z,Xs)  ■ f (Z,Xs-,x)dZdxs- 

1(6) 

In  the  above  equations,  Cj{  and  Yr  are  the  reaction  rates  and  species  mass 
fractions  obtained  from  the  steady-state  laminar  flamelet  calculations5.  The 
input  to  the  flamelet  calculations  includes  the  free-stream  values  of  species 
concentrations  (Yii,!#)  and  temperatures  (Ti,T2),  and  the  stoichiometric 
dissipation  rates,  which  have  to  be  modeled  from  the  LES  calculations.  The 
normalized  parameters  Da,  Ze  and  Ce  are,  respectively,  the  Damkholer 
number,  Zeldovich  number,  and  the  Heat  Release  parameter. 

The  joint  probability  density  function  of  the  mixture  fraction  and  its  dis- 
sipation rate,  f (Z,  xs),  contains  the  statistical  information  on  the  flamelets 
in  a turbulent  flow.  Statistical  independence  is  assumed  for  the  mixture 
fraction  (Z)  and  its  scalar  dissipation  rate  (x«): 

f (Z,  Xs ; x)  = f(Z;  x)  • /(x«;  x). 


As  most  reactions  occur  around  the  flame  sheet,  which  is  close  to  the  stoi- 
chiometric surface  in  a statistical  sense,  it  is  reasonable  to  assume  that 

/(x«;x)  = s (xst),  (2) 


and  /(Z;x)  is  usually  modeled  by  the  /3— form  PDF: 

ga—  1 /I  ry\b—  1 


f(Z-,x)  = 


(1  - Z)° 


B (a,  b) 


(3) 


where  a — Z ■ 


IM 1 

<z"2> 


1 


b=  a(l-  1),  and  B [a,  b)  = rSmih-  This 


r(o)+r(6)  • 

paper  compares  the  performance  of  (3— form  PDF  model  with  that  of  a 
Gaussian  PDF: 


f(Z;x) 


2tt  < Z”2  > 


\ 1/2 

t 

j exp 

— 

2 < Z"2  > 

/ 

(4) 


302 


XIAODAN  CAI  AND  FOLUSO  LADEINDE 


Similarity  assumption3 * * * 7  has  been  used  to  model  the  filtered  mixture 
fraction  variance  and  dissipation  rate  for  the  a posteriori  test  which,  re- 
spectively, are 


Z"2  = Z2-Z2~Ci 


and 


- ~ ndZdZ 

X — @2  ' D~ 


dxi  dx{ 


The  symbol  implies  test-level  filtering.  The  constants  C\  ~ 1.3  and 
C2  — IT  have  been  chosen  from  direct  numerical  simulations  that  were 
carried  out  as  part  of  the  present  study.  For  the  a priori  test,  they  are 
obtained  directly  from  the  DNS  data.  Cook  and  Riley7 *  proposed  the  use 
of  the  counterflow  assumption  to  close  the  conditional  filtered  dissipation 
rate  ( Xst ) with  the  above  filtered  dissipation  rate,  i.e., 


Xst  ^ X 


exp 

{—2  [er/c-1  (2Zst)]2j 

1 

So  exP  { 

-2[erfc~1(2Z)}2}f(Z)dZ 

(5) 


However,  it  is  argued  that  the  counterflow  structure  is  rarely  found  in  realis- 
tic turbulent  reacting  flows.  For  example,  recent  DNS  of  turbulent,  reacting 
mixing  layers22  found  that  the  chemical  reaction  occurs  typically  in  a shear- 
type  stretching  mode  instead  of  a counterflow  structure,  even  though  the 
concept  of  laminar  flamelets  was  still  applicable  in  this  case.  A mapping- 
closure  approach  has  been  attempted  by  the  current  authors23  to  replace 
the  counterflow  model.  However,  little  improvement  has  been  found.  One 
phenomenological  model  that  has  been  used  in  RANS  modeling  is  evalu- 
ated in  this  paper  but  within  the  framework  of  the  flamelet-LES  procedure 
being  reported  on.  The  phenomenological  model  takes  the  form 


Xst  - C3x. 


(6) 


3.  Numerical  Procedures 

Both  DNS  and  LES  are  performed  to  evaluate  the  subgrid  models  through  a 

priori  and  a posteriori  testings.  The  numerical  methods  employed  to  solve 

the  DNS/LES  equations  are  the  compact  schemes  for  the  spatial  derivatives 

and  the  classical  fourth-order  Runge-Kutta  scheme  for  the  time  integra- 

tion. This  is  combined  with  a high-order  filtering  procedure10  in  order  to 

suppress  numerical  instabilities  arising  from  the  unsolved  scales,  mesh  non- 

uniformities and  boundary  conditions.  It  must  be  noted  that  this  filtering 
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operation  has  a much  weaker  attenuation  effect  on  the  Fourier  amplitudes 
of  (pi  than  the  LES  filters  so  that  this  effect  won’t  mask  the  the  effect  of  the 
LES  filter.  The  LES  filter  scheme  uses  the  Simpson’s  averaging  scheme.  To 
account  for  the  boundary  effects,  the  Navier-Stokes  Characteristic  Bound- 
ary Condition  (NSCBC)14  is  extended  to  the  generalized  curvilinear  coor- 
dinate system  in  the  present  work.  The  steady  flamelet  equations  are  solved 
by  the  Newton- Raphson  method,  as  described  in  Ladeinde  et  allb. 

A single-step,  irreversible  chemical  reaction  of  the  type 

F + rO  — > (r  + 1)  P,  (7) 

is  used,  where  r represents  the  mass  stoichiometric  ratio  of  oxidizer  to 
fuel,  and  is  related  to  the  stochiometric  value  of  mixture-fraction  by  Zst  = 
[1  + r ■ Ffp/Vo]-1-  For  methane/air  combustion,  which  is  represented  by 
Case  2,  r = 4 and  Zst  — 0.055.  The  reaction  rate  assumes  the  form 

ur  = p2  ■ Da  ■ Y\  Y2  ■ exp 

The  two-dimensional  calculations  in  this  work  are  able  to  capture  the  essen- 
tial features  of  the  reaction  mechanism.  The  test  cases  are  listed  in  Table 

1.  Case  1 and  Case  2 involve  decaying  homogeneous  turbulence  with  ini- 
tial Reynolds  number  of  Re 0 = 250.  The  initial  field  for  Z is  random  with 
a pseudo  double-delta  probability  distribution17.  The  initial  fields  for  the 
mass  fractions  of  the  reacting  species  are  then  generated  from  Z assuming 
the  fast-limit  reaction.  The  smallest  turbulent  scales  under  the  Reynolds 
number  of  250  are  fully  resolved  by  the  DNS  grids  for  Case  1 and  Case  2. 
Case  3 involves  a spatially-developing  mixing  layer  with  an  inflow  Reynolds 
number  of  720  (based  on  vorticity  thickness).  The  convective  Mach  number 
is  0.125  with  a ratio  of  slow-  to  fast-stream  velocity  0.5.  Both  the  initial 
velocity  field  and  the  initial  mixture-fraction  field  use  the  hyperbolic  tan- 
gent profile,  and  the  initial  mass  fraction  fields  for  the  reacting  species  are 
generated  from  Z,  assuming  the  fast-limit  reaction  distribution.  The  inflow 
conditions  are  generated  by  superposition  of  small  perturbations  on  the 
mean  field.  To  facilitate  the  formation  of  roll-up  structures,  the  perturba- 
tions are  generated  from  the  fundamental  modes  (uq  = 1.3198)  obtained 
from  linear-stability  analysis24.  The  DNS  grids  for  Case  3 have  a resolution 
of  25ry  {rj  is  the  Kolmogorov  length  scale),  which  seems  fine  enough  for  the 
mixing  layer  turbulent  flows20  when  the  large-scale  structures  dominate  the 
flow  behaviors.  LES  was  carried  out  for  Case  1 and  Case  3 but  not  for  Case 

2.  However,  Case  2 was  used  in  the  evaluation  of  the  flamelet  model  (Figure 
!)• 
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Table  1.  Test  Cases. 


Test 

Flow  Type 

Da 

Ze 

Ce 

DNS  Grid 

LES  Grid 

Case  1 

Homogeneous 

0.5 

10 

0 

0 

133x133 

35x35 

Case  2 

Homogeneous 

0.055 

10 

0 

0 

133x133 

N/A 

Case  3 

Mixing  Layer 

0.5 

10 

5 

1.0 

375x99 

186x49 

4.  Evaluation  of  Models 

Flamelet  Calculation:  Prior  to  the  a priori  and  a posteriori  analyses  of 
model  performance,  flamelet  calculations  are  carried  out  to  obtain  the  val- 
ues of  u>i  and  F*  for  the  integrals  in  Eq.  (1).  The  calculations  are  performed 
in  the  mixture-fraction  space15.  Note  that  the  parameters  Da,  Ze  and  Ce, 
are  required  for  the  flamelet  calculations.  The  values  of  oij  and  Yt  from  the 
flamelet  calculation  are  then  constructed  in  the  form  of  a look-up  table  (in- 
dexed by  Z,  < Z"2  >,x  ) to  facilitate  the  LES  calculations.  Figure  1 shows 
the  product  mass  fraction  results  from  the  flamelet  calculations  for  Case 
1 and  Case  2.  The  effects  of  conditional  scalar  dissipation  rate  on  species 
distribution  appear  to  be  significant  for  the  chemistry  model  in  Eq.  (7), 
especially  around  the  stoichiometric  values  of  the  mixture-fraction.  This 
result  implies  that  the  scalar  dissipation  rate  model  has  significant  effects 
on  the  filtered  reaction  rates  and  filtered  mass  fractions  of  product  around 
their  maximum  values,  which  occur  close  to  the  stoichiometric  value  in 
mixture  fraction  space.  However,  it  is  noted  that  the  same  chemistry  model 
has  been  used  by  DesJardin  and  Frankel6  but  with  the  erroneous  assump- 
tion that  the  concentrations  of  the  species  were  independent  of  the  scalar 
dissipation  rate. 

a priori  analysis:  Data  sets  from  the  DNS  for  the  three  cases  are  used  to 
evaluate  the  accuracy  of  the  flamelet  models.  To  proceed  with  the  analysis, 
the  DNS  data  is  filtered  by  the  Simpson’s  scheme  onto  coarse  grid-points 
(A dns  — C ■ Ac,  where  A dns  is  DNS  grid-spacing  and  Ac  is  the  coarse- 
grid  spacing).  The  “exact”  values  for  Yp,  Co/ , Z,  Z”2  and  x are  used  to 
represent  their  corresponding  filtered  quantities  on  the  coarse  grid  point. 
The  last  three  quantities  are  then  used  to  obtain  the  model  values  for  Yp 
and  uj f using  the  flamelet  model  described  above.  Since  the  grid-coarsening 
factor,  C,  determines  the  number  of  sampled  DNS  points  in  each  LES  grid 
point,  it  may  affect  the  statistical  properties  of  the  embedded  flamelets. 
Therefore,  we  test  two  levels  of  this  parameter:  C — 2 and  C = 4. 

Figure  2 shows  the  contour  maps  of  product  mass-fraction  obtained 
from  DNS  for  Case  1 (eddy-turn-over  time  of  4.0,  Figure  2(a))  and  Case  3 
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(evolution  time  of  208,  Figure  2(b)).  The  contour  maps  evidently  show  that 
the  fine  flow-structures  are  suitably  resolved,  and  therefore  can  be  used  to 
evaluate  the  flamelet  chemistry  model. 

It  is  noted  that  the  present  flamelet-LES  approach  involves  a series 
of  models,  among  which  are  the  models  for  the  PDF  of  mixture-fraction 
{f  (z))  and  the  conditional  filtered  dissipation  rate  ( Xst )•  We  test  three 
models  for  / (z):  a)  5— PDF  (/  (z)  = 4 (z),  b)  /3-PDF  (Eq.  3),  c)  Gaussian 
PDF  (Eq.  4);  and  two  models  for  Xst  (Eq.  5 and  Eq.  6).  The  errors  from 
these  models  have  been  evaluated  by  two  criteria,  and 

V (x model -X eXact)  ^ where  X represents  Y„  orw/.  As  Criterium  2 puts  more 

weight  on  the  large  values,  it  is  considered  to  be  more  suitable  for  evaluating 
the  model  performance  in  the  sense  that  the  model  has  more  significant 
effects  on  the  behaviors  of  Yp  and  ujj  around  their  maximum  values,  which 
are  discussed  at  the  beginning  of  this  section.  Figures  3(a),  3(b),  and  3(c) 
present  the  model  errors  for  up  for  Cases  1,  2 and  3,  respectively,  with  grid- 
coarsening factor  (7  = 2.  Figure  3(d)  shows  the  model  errors  for  u>f  with 
(7  = 4 (Case  1).  The  values  of  C3  for  the  counterflow  model  in  the  figures 
represent  an  averaged  quantity  over  the  whole  domain  since  the  counterflow 
model  for  Xst  produces  a local  C3  value  at  each  grid  point.  From  these 
figures  we  see  that  the  errors  for  the  Gaussian  PDF  model  and  the  (3— PDF 
are  much  smaller  than  those  for  4 -PDF.  The  improvement  in  the  Gaussian 
PDF  and  (3- PDF  prediction  increases  with  increasing  values  of  (7.  Since 
the  use  of  4-PDF  implies  that  there  is  no  local  PDF  model  for  the  mixture 
fraction  on  each  LES  grid  point,  the  improvement  in  the  Gaussian  PDF  and 
(3— PDF  predictions  suggests  that  the  use  of  some  subgrid  PDF  models,  even 
though  simple,  will  be  very  helpful  in  LES  calculations,  a point  consistent 
with  the  results  obtained  by  Jimenez  et  al .21  It  is  also  observed  that  the 
Gaussian  PDF  model  consistently  performs  almost  as  well  as  the  (3— PDF, 
although  the  latter  is  used  more  often  in  the  literature.  The  results  here 
suggest  that  the  use  of  the  Gaussian  PDF  is  equally  acceptable.  This  point 
is  important  since  the  Gaussian  function  is  much  cheaper  to  calculate  than 
the  (3  function.  Furthermore,  it  can  be  seen  from  Figure  3 that  the  Xst  model 
also  has  a strong  effect  in  comparison  to  the  Gaussian  and  (3— PDF  models. 
The  graphs  clearly  show  that  an  optimal  C3  value  exists  for  each  case. 
It  may  be  said  that  the  performance  of  the  counterflow  model  is  modest 
in  comparison  to  the  performance  of  the  phenomenological  models  in  Eq. 
(6).  This  point  is  strengthened  by  Figure  4,  which  shows  the  correlation 
between  the  exact  and  modeled  values  of  d>f  in  the  form  of  scatter  plots 
for  the  three  cases.  Here,  the  results  from  the  counterflow  model  for  Xst  are 
compared  with  those  from  Eq.  (6)  using  the  optimal  value  of  coefficient  <73. 

a posteriori  analysis:  The  a posteriori  tests  are  considered  to  be  the 
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ultimate  tests  of  model  performance  in  the  sense  that  the  model  is  evaluated 
only  after  it  has  been  implemented  in  the  large-eddy  simulation.  Figure  5(a) 
shows  the  correlation  between  the  DNS  and  LES  values  of  Yp  (Case  1)  in 
the  form  of  scatter  plots.  It  appears  from  this  figure  that  the  LES-flamelet 
procedure,  along  with  the  assumptions  in  Eq.  (3)  and  Eq.  (6),  is  capable  of 
generating  accurate  results.  Further  evidence  comes  from  a comparison  of 
the  growth  rates  of  momentum  thickness  Q in  the  DNS  and  the  LES-flamelet 
calculations  for  the  spatially-developing  mixing  layers  (Figure  5(b)).  The 
LES-flamelet  model  gives  better  results  compared  to  those  without  the 
model.  However,  owing  to  the  combined  effects  of  numerical  discretization 
in  space  and  time,  and  averaging13,  a posteriori  tests  of  the  mixing  layer 
flows  do  show  some  scatter  for  the  correlation  between  the  DNS  and  LES 
values  of  Yp. 


5.  Conclusions 


We  have  shown  satisfatory  performance  of  a flamelet-LES  procedure  for  de- 
caying homogeneous  turbulence  and  the  spatially-developing  mixing  layer. 
The  chemical  reaction  used  in  the  tests  has  a moderate  Damkholer  num- 
ber and  a reasonable  Zeldovich  number  but  with  weak  heat  release.  The 
laminar-flame  structures  of  the  chemical  reactions  are  strongly  dependent 
on  the  conditional  scalar  dissipation-rate  and  are  thus  very  sensitive  to  the 
model  for  this  quantity. 

We  use  the  a priori  analysis  to  compare  the  model  values  of  Yv  and 
iof  with  their  corresponding  “exact”  values  from  DNS  data  sets.  The  error 
analyses  prove  that  the  use  of  some  subgrid  models  for  the  PDF  of  mixture- 
fraction,  even  though  simple,  yields  accurate  LES  calculations,  which  is 
consistent  with  the  observations  of  Jimenez  et  al.21  We  also  see  in  this  study 
that  the  Gaussian  PDF  model  consistently  shows  a comparable  performance 
to  the  more  common  (3— PDF,  and  may  provide  an  alternative  approach. 
This  point  is  important  because  the  Gaussian  function  is  much  cheaper  to 
calculate  than  the  (3  function.  Furthermore,  it  may  be  said  from  the  error 
analysis  that  the  counterflow  model  shows  only  a moderate  performance 

The  a posteriori  analysis  confirms  the  observation  that  the  performance 
of  the  flamelet-LES  procedure  consistent  with  the  model  assumptions  in  this 
paper  appears  to  be  satisfactory.  However,  owing  to  the  combined  effects  of 
numerical  discretization  in  space  and  time,  and  averaging,  a posteriori  tests 
of  the  mixing  layer  flows  show  some  scatter  for  the  correlation  between  the 
DNS  and  LES  values  of  Yp. 
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Figure  1.  Flamelet  results  for  Case  1 (a)  and  Case  2 (b),  showing  the  effects  of  different 
values  of  Xst 


Homogeneous  Turbulence 


Mixing-Layer 


Figure  2.  Contour  maps  of  product  mass-fraction  at  time  of:  a)  t = 4.0  (Case  1);  b) 
t = 208  (Case  3). 
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Case  1 Case  2 


Figure  3.  The  effect  of  Coefficient  Cs  on  the  model  errors  of  w/  for  the  three  cases 
((a)-(c.)),  with  C — 2.  The  error  plot  in  Figure  3(d)  is  for  Case  1 with  C — 4. 


Figure  5.  a posteriori  tests  of  model  performance:  a)  the  correlation  between  DNS  and 
LES  values  of  Yv  in  the  form  of  scatter  plots  for  Case  I;  b)  comparison  of  the  growth 
rates  of  the  momentum  thickness  9 from  DNS  and  flamelet-LES  calculations. 


